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Abstract 



The general 4D rotation matrix is specialised to the general 3D rotation matrix by equating its aoo element 
to 1. Its associate matrix of products of the left-hand and right-hand quaternion components is specialised 
correspondingly. Inequalities involving the angles through which the coordinate axes in 3D space are displaced 
are used to prove that the left-hand and the right-hand quaternions are each other's inverses, thus proving 
the Euler— Rodrigues formula. 
' A general procedure to determine the Euler parameters of a given 3D rotation matrix is sketched. 

By putting aoo = — 1 instead of +1 in the general 4D rotation matrix one proves the counterpart of the 
Euler— Rodrigues formula for 3D rotoreflections. 

Keywords: Euler-Rodrigues formula, Euler parameters, quaternions, four-dimensional rotations, three- 
dimensional rotations, rotoreflections 
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10 Notations 



a, b, c, d In Section [3J Euler parameters of a 3D rotation 

A Arbitrary 4D rotation matrix 

q . M Associate of matrix A (defined in [MEBI 2005]) 

Ml, Mr Matrices representing left- and right-multiplication by a unit quaternion, respectively 

^3 . L = a + bi + cj + dk Unit quaternion appearing as a leftmost factor in quaternion multiplication 

j3 1 P = u + xi + yj + zk Arbitrary 4D point represented as a quaternion 

G . R = p + qi + rj + sk Unit quaternion appearing as a rightmost factor in quaternion multiplication 

' OXYZ 3D Cartesian coordinate system 

• i-h ■ p Distance of point in 3D space from Z axis 

X' 
h : 

.5^,1 The quaternion representation theorem for 4D rotations 

Theorem: Each 4D rotation matrix can be decomposed in two ways into a matrix representing left -multiplication 
by a unit quaternion and a matrix representing right-multiplication by a unit quaternion. These decompositions 
differ only in the signs of the component matrices. 

In quaternion representation: P' = LPR; in matrix representation: P' = MlMrP. 



Outline of proof: Let Ml, Mr (Eq. [T]) be matrices representing left- and right-multiplication by a unit 
quaternion, respectively. Then their product A = MlMr (Eq. [2]) is a 4D rotation matrix. 

Matrix Ml is determined by four reals a, b, c, d satisfying the relation a 2 + b 2 + c 2 + d 2 = 1. Likewise, matrix Mr 
is determined by four reals p, q, r, s satisfying the relation p 2 + q 2 + r 2 + s 2 = 1. 

The 16 products ap, aq, ar, as, . . . , dp, dq, dr, ds are arranged into a matrix M (Eq. |4j, which has rank one 
and norm unity (when considered as a 16D vector), and is easily expressed in the elements of A. In [MEBI 2005] 
matrix M is denoted as the associate matrix of A. 
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Conversely, given an arbitrary 4D rotation matrix A (Eq. [3]), one calculates its associate matrix M according to 
Eq. [5] in the hope that it is a matrix of products ap, aq, ar, as, . . . , dp, dq, dr, ds which are not all zero. This 
hope is vindicated by proving that M has rank one whenever A is a 4D rotation matrix. For this proof one needs 
the general theorem that complementary subdeterminants of a rotation matrix of any number of dimensions are 
equal. 

The proof is completed by observing that the sum of the squares of the elements of M is unity, and concluding 
that two pairs of quadruples of reals a, 6, c, d; p, q, r, s exist satisfying a 2 + b 2 + c 2 + d 2 = 1, p 2 + q 2 + r 2 + s 2 = 1 
and differing only in sign. 
A complete proof is given in [MEBI 2005] . 
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2 3D Rotations 

Any 3D isometry with a fixed point O has an invariant plane and an invariant axis which intersect at right angles 
at O. 

An orientation-preserving isometry with fixed point O is a 3D rotation through a certain angle a around its 
invariant axis in its invariant plane. (Theorem of Euler) 

An orientation-reversing isometry with fixed point O is a 3D rotation through a certain angle a in its invariant 
plane combined with a reflection in that plane. Hence the name rotoreflection for such an isometry. 

In this section we prove the Euler-Rodrigues formula for 3D rotations. This formula arises from the gen- 
eral 4D rotation formula by setting p — a, q = —6, r = — c, s = —d in Eq. [2] Expressing this in quater- 
nion terms one obtains the Hamilton-Cayley formula P' = QPQ^ 1 , or expanded: u' + x'i + y'j + z'k = 
(a + bi + cj + dk) [u + xi + yj + zk) (a — bi — cj — dk). 

Proof: The general 3D rotation matrix can be obtained from the general 4D rotation matrix by putting aoo = 1, 
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oio = a 2Q — a 3 o = 0, agi = ap2 = ao3 = 0. Doing so one obtains 
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Its associate matrix reads 
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Comparing the off-diagonal elements with the corresponding elements in Eq. U one infers 

aq = —bp, ar = —cp, as = —dp, 
cq = br, dq = bs, dr = cs. 
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Let a be the rotation angle of the 3D rotation. The trace of its matrix is an + 022 + 033 = 2 cos a + 1, and the 
diagonal elements of the assoicate matrix are 



ap = (on + o 2 2 + 033 + l)/4 = (2cosa + 2)/4 



(10) 



bq = (an + a 2 2 + a 33 - 1 - 2an)/4 
cr = (an + a 2 2 + a 33 - 1 - 2a 2 2)/4 
ds = (an + a 22 + a 33 - 1 - 2a 33 )/4 



(2 cos a - 2au)/4 
(2 cos a - 2a 22 )/4 
(2 cos a - 2a 33 )/4 



(11) 
(12) 

(13) 



From Eq. [10] one deduces 



ap > 0. 



(14) 



The quantities an, 022, a 33 are the cosines of the angles through which the +OX, +OY and +OZ half-axes are 
displaced. These angles are at most a, therefore their cosines are at least cos a. For a proof of this see Section [5] 
We conclude that 

bq < 0, cr < 0, ds < 0. 



We recall the conditions 



From Eq. [8] follows the equality 



b 2 + c 2 +d 2 = l, p 2 



2,2,2 
7 +r + s 



a 2 (q 2 +r 2 + s 2 )=p 2 (b 2 +c 2 +d 2 ), 



which is because of Eq. [T|5]the same as a 2 (l — p 2 ) = p 2 (l — a 2 ), or equivalently: a 2 = p 2 . 
Taking into account Eq. [H] we conclude a = p. 

Substitute a for p in Eq. [8] to obtain aq = —ab, ar = —ac, as = —ad. 

If a ^ then q = —b, r — — c, s = —d. 

If a = then also p = 0. In that case one has 



c 2 + d 2 



1, q 2 + r 2 +s 2 



1. 



(15) 

(16) 
(17) 



(18) 



Eq. [5] is equivalent to (6,c, d) oc (q, r, s). Combining this with Eq. IT51 and Eq. [15] we conclude that q = —b, 
r = — c, s = —d also for a = 0. End of proof. 



3 



3 The Euler— Rodrigues formula 



By substituting a, —b, — c, — d for p, q, r, s in Eq. [2] and deleting the leftmost column and uppermost row one 
obtains the Euler-Rodrigues formula in its usual form: 



with the relation a 2 



a 2 + b 2 -c 2 - d 2 
2ad + 2bc 
-2ac+2bd 
b 2 + c 2 + d 2 = l. 



-2ad + 2bc 
a 2 -b 2 + c 2 - d 2 
2ab + 2cd 



2ac + 2bd 
-2ab + 2cd 



(19) 



This formula was first discovered by Euler ([EULE 1770]) and later rediscovered independently by Rodrigues 
([RODR 1840]), hence the terms Euler parameters for the parameters a, b, c, d and Euler- Rodrigues for- 
mula (or formulae) for the rotation matrix in terms of a, 6, c, d. 



3.1 Determination of the Euler parameters for a given 3D rotation 

By substituting a, —b, — c, — d for p, q, r, s in the associate matrix (Eq. or directly from the Euler- 

Rodrigues formula and the accompanying relation, one obtains the ten equations below in the four unknowns 
a, b, c, d: 

a 2 = (l + an+a 22 + a 33 )/4, (20) 

b 2 = (1 + aii -a 22 -o 33 )/4, (21) 

c 2 = (1 - an + a 22 - a 33 )/4, (22) 

d 2 = (1 - ou - 02a + oas)/4, (23) 



ab = (a 32 - a 23 )/4, (24) 

ac = (ois - a 3 i)/4, (25) 

ad = (oai-oi 2 )/4, (26) 

cd = (a 32 +a 23 )/4, (27) 

db = (oi3 + o 3 i)/4, (28) 

be = (a 2 i+ai 2 )/4, (29) 



which are consistent and admit a pair of real solutions if and only if the given matrix is indeed a 3D rotation matrix. 

Assuming that the given matrix an . . . a 33 is indeed a 3D rotation matrix, one first determines the absolute 
values of a, 6, c, d from Eqs. [201 ■ ■ ■ 1231 next one arbitrarily fixes the sign of one of the nonzero quantities a, 6, 
c, d (at least one of them is nonzero) ; finally one determines the signs of the other ones from the remaining six 
equations. 

There always exists a subset of Eqs. [20] . . . [29] from which one can determine a, b, c, d save for their signs; the 
remaining equations provide opportunities for cross-checking the calculations and for averaging numerical errors. 
Think in this connection of ill-conditioned data and of input rotation matrices affected with errors. 
In practice one uses Eqs. [M] . . . [H] for rotations through small angles, and Eqs. [23 . . . [2H1 for rotations through 
angles close to it. 

These procedures are well-known in the fields of aviation and astrodynamics; see for instance [BATT 1999] and 
[SHEP 1978]. 



4 3D Rotoreflections 

In this section we prove the counterpart of the Euler- Rodrigues formula for 3D rotoreflections. This formula 
arises from the general 4D rotation formula by setting p = —a, q = b, r = c, s = d in Eq. [5] Express- 
ing this in quaternion terms one obtains the Hamilton— Cayley— like formula P' = —QPQ , or expanded: 
u' + x'i + y'j + z'k = — [a + hi + cj + dk)(u + xi + yj + zk)(a — bi — cj — dk). 
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Proof: The general 3D rotoreflection matrix can be obtained from the general 4D rotation matrix by putting 
aoo = — 1, »io = a 20 = a 30 = 0, aoi = ao2 = 003 = 0. Doing so one obtains 
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Comparing the off-diagonal elements with the corresponding elements in Eq. [4] one infers 

aq = —bp, ar = —cp, as = —dp, 
cq — br, dq = bs, dr = cs. 



(32) 
(33) 



Let a be the rotation angle of the 3D rotoreflection. The trace of its matrix is an + 022 + a 33 = 2 cos a — 1, and 
the diagonal elements of the assoicate matrix are 



bq 

cr 
ds 



ap = (on + a22 + 033 - l)/4 = (2cosa-2)/4 



(on + a 22 + a 33 + 1 - 2a n )/4 
(an + a 22 + a 33 + 1 - 2a 22 )/4 
(an + a 22 + a 33 + 1 - 2a 33 )/4 



(2 cos a - 2a n )/4 
(2 cos a - 2a 22 )/4 
(2 cos a - 2a 33 )/4 



(34) 

(35) 
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From Eq. [34] one deduces 



ap < 0. (38) 



The quantities an, 022, a 33 are the cosines of the angles through which the +OX, +OY and +OZ half-axes are 
displaced. These angles are at least a, therefore their cosines are at most cos a. For a proof of this see Sectional 
We conclude that 

bq > 0, cr > 0, ds > 0. (39) 

We recall the conditions 

,2 , u2 , „2 , d 2 = ^ p 2 + ? 2 + r 2 + s 2 = j ( 4Q ) 



a + + c 



From Eq. [32] follows the equality 

a 2 {q 2 +r 2 + s 2 )=p 2 {b 2 +c 2 +d 2 ) 1 (41) 

which is because of Eq. [40] the same as a 2 (l — p 2 ) = p 2 (l — a 2 ), or equivalently: a 2 = p 2 . 

Taking into account Eq. [38] we conclude a = —p. 

Substitute —a for p in Eq. [32] to obtain aq = ab, ar = ac, as = ad. 

If a 7^ then q = b, r = c, s = d. 

If a = then also p = 0. In that case one has 

b 2 + c 2 +d 2 = l, q 2 + r 2 + s 2 = 1. (42) 

Eq. [33] is equivalent to (6, c, d) oc (q, r, s). Combining this with Eq. 1421 and Eq. [39] we conclude that q = b, r = c, 
s = d also for a = 0. End of proof. 
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The counterpart of the Euler-Rodrigues formula for 3D rotoreflections reads 



-a 2 -b 2 + c 2 +d 2 
-2ad - 2bc 
+2ac - 2bd 

with the relation a 2 + b 2 + c 2 + d 2 = 1. 



+2ad - 2bc -2ac - 2bd 

-a 2 + b 2 -c 2 +d 2 +2ab - 2cd 

-2ab - 2cd -a 2 + b 2 + c 2 - d 2 



(43) 



5 Proofs of the angle inequalities for 3D rotations and rotoreflections 



The matrix of a 3D rotation around the Z axis reads 

cos a —sin a 
R = sin a cos a 
1 

Let (j be the angle through which the ray from O through a point (x, y, z) is displaced. An easy calculation yields 
(putting p 2 = x 2 + y 2 ) 

cos (3 = (p 2 cos a + z 2 )/(p 2 + z 2 ), 
cos/3 - cosa = (p 2 cosa + z 2 - p 2 cosa - z 2 cosa)/(p 2 + z 2 ) = z 2 (l - cosa)/(p 2 + z 2 ) > 0. 



The matrix of a 3D rotoreflcction with the Z axis as its axis reads 



R 



cos a —sin a 
sin a cos a 




Let (3 be the angle through which the ray from O through a point (x, y, z) is displaced. An easy calculation yields 
(putting p 2 = x 2 + y 2 ) 

cos/3 = (p 2 cos a — z 2 )/(p 2 + z 2 ), 
cos/3 — cos a = (p 2 cos a — z 2 — p 2 cos a — z 2 cosa)/(p 2 + z 2 ) = —z 2 (l + cosa)/(p 2 + z 2 ) < 0. 
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